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ABSTRACT 

SAX J2 103. 5 +4545 is the Be/X-ray binary with the shortest orbital period. It shows ex- 
tended bright and faint X-ray states that last for a few hundred days. The main objective of 
this work is to investigate the relationship between the X-ray and optical variability and to 
characterise the spectral and timing properties of the bright and faint states. We have found a 
correlation between the spectral and temporal parameters that fit the energy and power spec- 
tra. Softer energy spectra correspond to softer power spectra. That is to say, when the energy 
spectrum is soft the power at high frequencies is suppressed. We also present the results of 
our monitoring of the Ha line of the optical counterpart since its discovery in 2003. There is a 
correlation between the strength and shape of the Ha line, originated in the circumstellar en- 
velope of the massive companion and the X-ray emission from the vicinity of the neutron star. 
Ha emission, indicative of an equatorial disc around the B-type star, is detected whenever the 
source is bright in X-rays. When the disc is absent, the X-ray emission decreases significantly. 
The long-term variability of SAX J2103. 5+4545 is characterised by fast episodes of disc loss 
and subsequent reformation. The time scales for the loss and reformation of the disc (about 2 
years) are the fastest among Be/X-ray binaries. 
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1 INTRODUCTION 

High-mass X-ray binaries (HMXB) typically consist of a neutron 
star orbiting around an early-type star. They divide into two major 
groups according to the luminosity class of the optical component: 
supergiant X-ray binaries if the massive companion is an evolved 
star (luminosity class I-II) and Be/X-ray binaries (BeX) if it is a 
main-sequence or giant star. The vast majority of BeX are tran- 
sient systems, althoug h there is a small group of persistent sources 
jReig&Rochei ri999). Transient BeX have orbital periods in the 
range 20-100 days and spin periods shorter than ~200 s. In con- 
trast, persistent BeX present orbital periods typically longer than 
~100 days and spin periods longer than ~200 s. Supergiant X-ray 
binary systems have narrower orbits with orbital periods shorter 
than 10 days and contain slow-rotating neutron stars, with spin pe- 
riods of the order of minutes to hours. 

SAX 12103.5+4545 i s an u nusual HMXB. Its spin, P spm = 
358.6 s dHulleman et al.1 Il998h . and orbital, P orb = 12.7 d 
(Bavk aTet al.ll2002i) periods are typical of supergiant X-ray bina- 
ries. In fact, SAX 12103.5+4545 falls in the wind-fed supergiant 
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region dReig et alj2004l) of the P orb - P S pm diagram dCorbetlll986h . 
However, the primary c omponent of the b inary is a main-sequence 
star of spectral type B0 dReig et alj|2004l) . i.e. SAX 12103.5+4545 
belongs to the BeX category. SAX J2103.5+4545 alternates X- 
ray bright states that last for a few months with exte nded (as 
long as few years) X-ray faint states (Bavkal et al. 2002). During 
the bright states SAX J2 103. 5 +4545 shows large spin-up episodes 
dlnam et alj|2004l : Isidoli et al.ll2005l) . in dicating that an acc retion 
disc is formed around the neutron star dBavkal et al.ll2002T) . The 
scarce optical data seem to suggest that during the bright state the 
Ho- line turns into emission, implying that a decretion circumstellar 
disc is formed around the Be star dReig et alj|2004h . These discs 
are probably sh ort lived and appe ar during the high X-ray emis- 
sion states only. iReiget all d2005l) found that SAX J2103.5+4545 
was emitting X-rays even after the complete loss of the circumstel- 
lar disc. They argue that in this state the X-rays are the result of 
wind-fed accretion. 

Since the discovery of SAX J2103.5+4545 in February 1997, 
the system has been observed in the X-ray band on numerous oc- 
casions. The first X-ra y observations were made with BeppoSAX 
dHulleman et"ai]|l998h . The source was active for about eight 
months and reached a peak intensity of 20 mCrab (2-25 keV) on 
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Table 1. Ho- equivalent width measurements (Icr errors). See Fig.[3]for a 
representative example of the line profiles. 



Date Julian date EW(Hff) Ha Telescope 

(2,400,000+) (A) profile" 



a: ABS: absorption, CQE: absorption with central quasi-emission peak 
DPE: double-peak emission, SH: shell profile. 



April 11, 1997. The X-ray emission showed pulsations with a pe- 
riod of 358.61 s. 

The All Sky Mon itor on board RXTE detected a second bright 
state two years later dBavkal et al J 1 20001) . reaching a peak inten- 
sity of 27 mCrab (2-12 keV) on October 28, 1999. The continu- 
ous monitoring of SAX J2103.5+4545 by RXTE and INTEGRAL 
since July 2002 allowed detailed pulse freq uency analysis and 
the determination of the o rbital parameters dBavkal et alj 120071 : 
ICamero-Arranz et alj|2007l) : the system has a moderately eccentric 
orbit withe = 0.401±0.018, an orbital period of 12.66528±0.00051 
days, and a semi-major axis of 80.8 ± 0.7 lt-s. 

In addition to BeppoSAX and RXTE, SAX J2 103. 5 +4545 has 
been observed with XMM-Newton and INTEGRAL in the bright 
state, hence allowing the characterisation of the energy spectrum 
from 0.1 to 150 keV. At low energies (< 4 keV), the X-ray energy 
spectrum presents a soft spectral component, consistent with black- 
body emission w ith kT = 1.9 keV and emitting radius of ~0.3 km 
dlnam et alj2004h . The broad-band (4-150 keV) energy spectrum is 
well fitted by a power-law (r = 1-1.5) with an exponential cutoff 
(£ cut = 8- 19 keV), a nd with a Kg fluorescence emission line (~6.4 
keV) from cool iron telav et alj2004l : lsidoli et alj2005h . 

The sp e ctral analysis of RXTE data performed by 
Bavk al et al ] d2002h covering the interval November 1999- 
August 2000 indicate that the spectrum of SAX J2103.5+4545 
is harder during bright states. This result was confirmed by the 
X-ray colour analysis performed by ICamero-Arranz et al.1 (2007) 
using six years worth of RXTE data. Transient qu asi-periodic 
oscilla tions around 0.044 Hz have been reported by llnam et all 
J2004I) . 

The excellent monitoring of SAX J2 103.5+4545 in the X-ray 
band contrasts with the scarcity of data at other wavelengths. The 
only published optical /near-IR observa tions of SAX J2103. 5+4545 
are those reported by Reig et al. ( 2004:)). They identified the opti- 
cal counterpart with a moderately reddened V=14.2 B0V star, at a 
distance of 6.5 kpc. In this paper we present the results of our mon- 
itoring of the Ha line for the 2003-2009 interval and investigate 
the X-ray variability of SAX J2103. 5+4545 in correlation with the 
optical data. X-rays result from accretion of matter from the op- 
tical companion's circumstellar disc onto the neutron star surface, 
hence providing information about the physical conditions in the 
vicinity the compact object. The Ha line is formed in the Be star's 
disc as reprocessed radiation from its photosphere. Since the disc 
constitutes the reservoir of matter available for accretion a corre- 
lated X-ray/optical study is the best way to unveil the nature of the 
variability in this system. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Optical spectra 

We have been observing SAX J2103. 5+4545 regularly since the 
identification of the optical counterpart in 2003. We have analysed 
all available optical o bservations including those already published 
bv lReigetal] d2004l) . The new optical spectroscopic observations 
were obtained mainly from two sites: the Skinakas (SKI) obser- 
vatory in Crete (Greece) and from the Fred Lawrence Whipple 
(FLW) observatory at Mt. Hopkins (Arizona). In addition, SAX 
J2103. 5+4545 was observed in service time from the William Her- 
schel Telescope (WHT) at the observatory of El Roque de los 
Muchachos in La Palma (Spain) on 1 September 2004. One more 
spectrum was made from the 2.2-m telescope at the observatory of 
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Table 2. Optical photometric observations of SAX J2103. 5+4545. 



Date Julian Date B V R I 

2,400,000+ mag mag mag mag 



08-06-2003 
24-08-2003 
20-05-2004 
05-07-2004 
24-08-2004 
14-09-2004 
01-10-2004 

26- 06-2005 

27- 07-2005 
20-08-2005 
18-08-2006 
16-07-2007 

01- 09-2007 

02- 09-2007 
26-10-2007 
05-08-2008 



52799.467 
52876.408 
53146.558 
53192.354 
53242.477 
53263.375 
53280.370 
53548.497 
53579.502 
53603.414 
53966.519 
54298.493 
54345.426 
54346.496 
54400.242 
54684.434 



15.35±0.03 
15.36±0.02 
15.41±0.01 
15.42±0.02 
15.39±0.02 
15.39±0.02 
15.41±0.02 
15.39±0.02 
15.38±0.02 
15.42±0.02 
15.42±0.03 
15.12±0.03 
15.21±0.03 
15.12±0.03 
15.37±0.03 
15.41±0.02 



14.22±0.02 
14.25±0.03 
14.32±0.01 
14.33+0.01 
14.31+0.01 
14.27+0.02 
14.31+0.02 
14.30+0.01 
14.29+0.01 
14.35+0.03 
14.32+0.03 
13.75+0.03 
13.94+0.02 
13.89+0.03 
14.21+0.03 
14.29+0.02 



13.48+0.02 
13.57+0.03 
13.63+0.01 
13.64+0.02 
13.61+0.01 
13.58+0.01 
13.61+0.02 
13.61+0.02 
13.60+0.01 
13.66+0.02 
13.63+0.02 
12.93+0.03 
13.15+0.03 
13.11+0.02 
13.57+0.04 
13.60+0.02 



12.85+0.03 
12.87+0.01 
12.92+0.02 
12.85+0.02 
12.84+0.03 
12.86+0.02 
12.86+0.04 



12.05+0.03 
12.32+0.04 
12.28+0.02 
12.78+0.05 
12.87+0.02 



Calar Alto (CAL) in Almerfa (Spain) on 3 December 2003. Tabled 
gives the log of the spectroscopic observations and the value of the 
Ho- equivalent width. 

The 1.3 m telescope of the Skinakas Observatory was 
equipped with a 2000x800 (15 /jm) pixel ISA SITe CCD and a 
1302 1 mnT 1 grating, giving a nominal dispersion of ~1 A/pixel. 
The instrumental set-up during the service WHT observation con- 
sisted of the ISIS spectrograph with the R1200R grating, giv- 
ing a dispersion of 0.23 A/pixel and covering the wavelength in- 
terval between 6040-6900 A. The spectrum was obtained with 
the 4096x2048 (13.5 /im) pixels MARCONI2 CCD. The reduc- 
tion of these spectra was made using the STARLINK Figaro and 
IRAF v2.14 packages, while their analysis was performed using 
the STARLINK Dipso package. The FLW observations of SAX 
J2103. 5+4545 were made in q ueue mode with the 1 .5-m telescope 
and the FAST-II spectrograph dFabricant et"al]f l998). Three obser- 
vations were obtained with the 1200 1 mm" 1 grating and one ob- 
servation (5 December 2007) with the 600 1 mm" 1 grating. The tilt 
of the grating was selected in order to cover the Ha line. The data 
were analysed with the RoadRunner package dTokarz & Rolll 19971) 
implemented in IRAF. Spectra of comparison lamps were taken be- 
fore each exposure in order to account for small variations of the 
wavelength calibration during the night. 

The field around SAX J2103. 5+4545 was observed through 
the Johnson B, V, R, and / filters from the 1.3 m telescope of 
the Skinakas observatory. Observations obtained before July 2007 
were performed using a 1024 x 1024 SITe CCD chip with a 24 
/jm pixel size (corresponding to 0.5" on sky). From July 2007 a 
2048x2048 ANDOR CCD with a 13.5 fxm pixel size was used 
(cor responding to 0.28" on sky). Standard stars from the Landolt 
list fcandoltll2009h were used for the photometric calibrations and 
transformation equations. Reduction of the data was carried out in 
the standard way using the IRAF tools for aperture photometry. 

2.2 X/y-ray observations 

We have analysed all archived observations of SAX J2103. 5+4545 
made by RXTEfPCA and INTEGRAL/IS GRI . The RXTEfPCA con- 
sists of five Proportional Counter Units (PCU), which are sensitive 



Table 3. Summary of the RXTE observations. 
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MJD 
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ID 
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time (ks) 


cs _1 
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51501.1846 


148.05 


20.1 
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24.3 


57 
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51606.2713 


154.16 


21.4 


79 
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51764.0044 


60.35 


7.0 


40 


P60409 


52019.0087 


5.20 
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41 


P70082 


52439.9567 


446.51 


33.2 


71 


P90097 


53029.1948 


252.86 


9.94 


67 


P92436 


54221.0628 


27.98 


40.54 


86 



a: 2-60 keV PCU2 background-subtracted. Time bin: 16 s 

b: root-mean-square: standard deviation over the mean intensity 



to X-rays in the 2-60 keV energy range dJahodaetal.ll 19961) . The 
INTEGRAL Soft Gamma Ray Imager (ISGRI) is the lower energy 
detector of the Imager on Board INTEGR AL Satellite (IBIS) and 
operates in the 15-1000 keV energy range dLebrum et alj2 003). 

In addition, data from two all sk y monitors, namely the 
All-Sky Monitor (ASM) onboard RXTE dLevine et al.l|l996T) and 
the B urst Alert Telescope (BAT) onboard SWIFT (Barthelmv 
2000) have been used to study the long-term variability of SAX 
J2103. 5+4545 (Fig. [TJ. These two instruments produce daily flux 
averages in the energy range 1.3-12.1 keV and 15-150 keV, respec- 
tively. 

RXTEfPCA data can be collected and telemetered to the 
ground in many different ways depending on the intensity of the 
source and the spectral and timing resolution desired. In this work 
we used the two standard modes: Standard! provides 0. 125-s reso- 
lution and no energy resolution; in the Standard! configuration data 
are accumulated every 16 seconds in 129 channels. For the tim- 
ing analysis GoodXenon data, providing maximum spectral (256 
channels) and temporal (2~ 20 s) resolution, have also been em- 
ployed. The number of active proportional counter units varied 
during the observations. However, PCU2 was always on. To avoid 
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Figure 1. Long-term evolution of the Ha equivalent width, R band and 
B — V colour. Also shown is the long-term variability in three X-ray bands, 
corresponding to the RXTE/ ASM, INTEGRALflSGSl and SWIFT/BAT light 
curves. 

complications with the different response of the each unit and the 
varying number of active detectors we used PCU2 data for the 
spectral analysis. Note that the power spectral analysis is not af- 
fected by the number of active units as we used the rm s normali- 
sation ( Bello ni & Hasingerll99d : lMivamoto et all 199 ll) , where the 
power is normalised to the mean intensity of short (128 s) segments. 
Table[3]gives a log of the RXTE observations. 

A total of 3513 public data pointings of INTEGRAL included 
SAX J2103.5+4545 in the IBIS/ISGRI field of view. However, 
only in 269 the source was detected with significance greater than 
3cr above the noise level (detection level greater than 7 based on 
the standard Offline Scientific Analysis (OSA) softwar^H), while in 
1436 additional observations we obtained marginal or non reliable 
detections. In the remaining 1808 pointings the source was below 
the sensitivity level of IBIS/ISGRI. 

OSA version 7.0 have been used to perform light curve and 
spectra extraction for each science window. To obtain average spec- 
tra corresponding to several science wind ows we used an altern a- 
tive event extraction method developed by Ferrigno e t alj d2007t) . 



3 RESULTS 

3.1 Optical photometry 

FigureQjshows the evolution of the Ho- equivalent width (EW(Ha)), 
R-band magnitude and the B - V colour together with the long- 
term X-ray variability in three different energy bands. Despite the 
observational gaps in the optical data a correspondence between the 
variability in the optical band and the X-ray intensity is apparent. 
Bright X-ray states correspond to enhanced optical activity, both in 

1 The OSA software generates sky images and searches for significant 
sources. A source is considered as detected if it has a detection level greater 
than 7, which is equivalent to the classical 3cr level above the measured 
noise. 
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Figure 2. EW(Hci) and I-band magnitude as a function of ASM intensity. 



the photometric magnitudes and colours and the strength of the Ha. 
This correlation is more clearly seen in Fig. [2] where the EW(Ho') 
and I-band magnit ude as a function of A SM intensity is displayed. 

As expected jDoughertv et aljfl994l and references therein), 
the contribution of the Be star's disc to the photometric magnitudes 
increases with wavelength. As can be deduced from Table [2] the 
amplitude of variability of the optical magnitudes between maxi- 
mum and minimum brightness is 0.3 mag in the B band, 0.6 mag 
in the V band, 0.7 mag i n the R band and 0.8 mag in the I band. 

In lReig etalj 1:2004) values of the reddening and distance were 
derived from the spectra and the photometric observations. Since 
i) a much longer set of optical photometric observations is now 
available and ii) a longer extended low optical state is observed, 
where presumably the underlying B star is exposed and photo- 
spheric emission without substantial contribution from the disc is 
detected, it is justified to revisit these values. The main source of 
uncertainty in estimating the distance stems from the uncertainty 
of the intrinsic colours and absolute magnitudes associated to each 
spectral and luminosity class. Difference between different cali- 
brations m ay amount up to 1.5 mag in the values of the absolute 
magnitude (Weg nerhoOrj) . while differences of 0.04 mag in intrin- 
sic (B — V)p are fou n d in studies from d i fferen t authors jjohnsonl 
1 1 9661 : IWegnerj[l994h . IJaschek & Gomezl Jl998h analysed the ab- 
solute magnitude of about one hundred MK standards and con- 
cluded that the intrinsic dispersion of the mean absolute magnitude 
amounts to 0.7 mag. We have adopted this value as the error on the 
absolute magnitude. 

To estimate the distance to SAX J2103. 5+4545 we use the ob- 
servations showing the bluest colours (those from 20 August 2005, 
MJD 53603), as no contribution from the circumstellar disc is ex- 
pected. By comparing the observed colour (B —V) = 1.0 7 ± 0.03 
with t he expected one of a BOV star (B - V) = -0.28 jjohnsorj 
1 19661 : [Gutierrez-Morendll979l : fWegner 1994), we derive a colour 
excess of E(B - V) = 1.35 + 0.03. Taking the standard reddening 
law A v = (3.1 ± 0. l)E(B - V) and assuming an average absolute 
magnitude for a BOV star of M v = -4.0 ± 0.7 dVacca et alll 19961: 
IWegnerll2006l) the distance to SAX J2103.5+4545 is estimated to 
be ~ 6.8±2.3 kpc. The final error was obtained by propagating the 
errors of (B - V) (0.03 amg), A v (0.17 mag) and M v (0.7 mag). The 
estimated distance is consistent with the value of 4.5+0.5 kpc de- 
rived from the rela tionship between th e pulse frequency derivative 
and the X-ray flux ( iBavkal et ai]|2007l) . 
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Figure 3. Ha line profiles at different X-ray activity states, black: emission (DPE), blue: absorption (ABS), orange: central quasi-emission (CQE) and cyan: 
shell profiles (SH). Ha- emission occurs when the source is X-ray active. See the electronic edition of the Journal for a colour version of this figure 



3.2 Ha line profiles 

The monitoring of SAX J2103. 5+4545 in the optical band reveals 
that the Ho- line is highly variable, both in strength and shape. In 
Be stars, the strength and shape of the Ho- line provide information 
about the physical conditions in the circumstellar disc. The mea- 
surements of the EW(Ha) are given in Table [T] A negative value 
indicates that the line appears in emission. The fourth column of 
Table Q] indicates the profile of the line. We have observed four 
different profiles (see Fig. [3J: emission with a shell profile (SH), 
absorption (ABS), emission with a split profile (DPE) and central 
(quasi)-emission peak (CQE). The difference between double-peak 
profile and a shell profile is simply that in the latter the depression 
between the two peaks extends below the stellar continuum. CQE 
is a type of absorption line profile in which t he central part of th e 
line core exhibits a weak local flux maximum teivinus etal.1 19991) . 

What determines whether the Ho- line appears in absorption 
or emission is the presence of an equatorially concentrated circum- 
stellar disc around the Be star. The absorption lines are formed in 
the photosphere of the star (i.e., no disc is present) while the emis- 
sion profile is the result of recombination radiation from ionised 
hydrogen in the hot, extended circumstellar envelope surrounding 
the central Be star. CQEs appear when the innermost regions of the 
disc are being supplied with matter. The deep central absorption 
in a shell profile is merely due to self-absorption in t he envelope 
as a c onsequence of a high inclination angle (see e.g. iHanuschikl 
1 19951) . CQE and shell features are related phenomena. CQEs oc- 
cur only when shell lines are also present and the lines showing 
CQE profiles are those also showing a shell feature. Figure[3]show 
a representative example of each one of these line profiles. 

During the bright states, the Ho- line exhibits an emission pro- 
file, indicating the presence of a circumstellar disc, whereas during 



the faint states the source shows a purely photospheric absorption 
profile or a CQE-type profile. 



3.3 X-ray states 

The long-term X-ray variability of SAX J2103. 5+4545 is charac- 
terised by bright and faint states (Fig.Q. The bright states appear as 
outbursts that last typically for 15-20 orbital cycles, i.e. a few hun- 
dred of days. These outbursts do not occur in a predictable manner 
but they do show a general repeatable pattern with regard to the 
duration and profile of the outburst. In general the outburst begins 
with a sharp flare that last for one or two orbital periods. This flare 
is followed by a progressive increase in the X-ray intensity until 
a maximum is reached at about one order of magnitude above the 
quiescent level. Assuming an absorbed power-law with exponential 
cutoff (see below) and a distance of 6.8 kpc, the X-ray luminosity 
at the peak of the outbursts ranges between 0.6 - 1.0 x 10 37 erg s" 1 , 
while at the peak of the flare the luminosity is typically a factor of 2 
higher. During the bright state a modulation of the X-ray intensity 
with orbital phase is clearly seen (Fig. |4j, but it is not detected in 
the faint state. The average intensity level of the faint state is ~ 5 c 
s _1 PCU', which corresponds to a X-ray luminosity of 3.3 x 10 35 
erg s _I in the 3-30 keV range. Although the count rate may go down 
to ~ 1 c s PCU~' in the faint state, the source is detected in all 
PC A pointings. 

In order to characterise the spectral and timing properties of 
the source at different states, we selected two subsets of PCA data 
corresponding to the bright state and two to the faint state. Our se- 
lection was based on the sour ce activity, i.e. its PCA flux, and pul- 
sar spin-up rate (see Fig. 12 in lCamero-Arranz et all2007t) . During 
the faint state the spin-up rate is close to zero (or even negative), 
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Figure 4. Bottom panel: RXTE/ASM light curve (1.3-12.1 keV) since April 1996. A bin size equal to the orbital period, P ol b = 12.67 d is used. Middle panel: 
RXTE/PCV2 (2-60 keV) light curve showing bright (Bl, B2) and faint states (Fl, F2). Top panel: Detailed view of the PCU2 light curve. The 12.7-d binned 
PCA light curves is superposed on the average count rate per observation interval (typically, a few thousand seconds integration time). Note the modulation of 
the X-ray intensity with the orbital period during the bright states and the lack of it during the faint state. 



whereas during the bright state it is >2 x 10~ 13 Hz s -1 . Thus, faint 
states correspond to epochs of little or no mass accretion rate. The 
duration (typically 8-10 orbital periods) of the states was chosen 
such that it includes the rising and declining phase of the bright 
state and provides good statistics for the faint state. The time ranges 
of these intervals are given in Table 4 in Modified Julian Date 
(MJD), below the state names. Additionally, we selected a one-day 
subset of PCA data during the flare at MJD 54221 (Fig.gJ. Note, 
however, that RXTE/PCA observations were triggered by the sud- 
den increase of the SAX J2103. 5+4545 flux in the ASM light curve. 
Therefore, the PCA observations did not catch the peak of this flare 
but began when the flux was already in the declining phase. For 
each one of these states an average energy spectrum was obtained. 
A systematic error of 0.6% was added in quadrature to the statisti- 
cal error in the average energy spectra. 

An absorbed power law with a high-energy cut off and an iron 
line at ~ 6.4 keV provided good fits to all energy spectra in the 
range 3-30 keV. The spectra of the bright state are distinctly harder 
than those of the faint state, in agreement with previous findings 
(Bavk aTet al.l2002llCamero-Arranz et al.l2007l) . The cutoff energy 
in the faint states occurs at slightly higher energies than in the bright 
state, whereas the central energy of the iron line is slightly larger 
in the faint state. The line width, however, is not well constrained. 



If let as a free parameter then the line is broader in the faint state. 
However, equally good fits are obtained in a relatively wide range 
of values (cr = 0.3 - 1.0 keV) if this parameter is fixed in the faint 
state. While the addition of an edge at ~9 keV clearly improves the 
fit in the bright stat^HL no such edge is required for the faint state. 
The value of the best-fit parameters of the flare spectrum agree with 
those of the bright state within the errors. 

We also obtained an INTEGRAL/IBIS average energy spec- 
trum for the bright2 and faint2 states covering the energy range 
20-100 keV. We selected science windows where a statistically sig- 
nificant detection of SAX J2103. 5+4545 was found and in which 
the source lied within 12° off the centre. This radius includes all 
observations in which the source was within the PCFOV but it dis- 
regards the observations in which the source was too close to the 
detector edge. Figure[5]shows the joint average PCA+ISGRI spec- 
tra corresponding to the flare, bright and faint states. To fit these 
spectra we used the same model as for the PCA data simply adding 
one extra cutoff power law and a relative normalisation factor (fixed 
to 1 for the PCA spectrum). However, we fixed the photon index of 



2 The probability that the improvement of the fit in the bright state occurs 
by chance by the addition of this component is ^10~ 3 based on the F-test. 
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Table 4. Spectral parameters resulting from the fits to the RXTE/PCA alone and RXTE/PCA+INTEGRAL/ISGRl energy spectra. Errors denote 90% confidence. 



Parameter 


Bright 1 
51559-51669 


Bright2 
52602-52810 


Faint 1 
51722-51880 


Faint2 
53180-53348 


Flare 
54221 






RXTE/PCA enen 


;y spectra 






N" 
H 

Photon index 
£cut (keV) 
PL norm* 
£ Fc (keV) 
o-Fe (keV) 
EW(Fe) (eV) 
Eedge (keV) 
r 

Flux c (3-30 keV) 


3.7±0.5 

0.85±0.05 

13 9 +09 
_ Q2 

2 30+ 09 
z - -0.04 

6.44±0.08 
0.5±0.2 
160±30 
9.7±0.4 
0.05±0.01 
4.8 
1.5/45 


3.7±0.4 
0.81±0.02 

13.7±0.3 

•j £+0.09 
J -0.06 

6.46±0.08 
0.4±0.1 
120±20 
9.7±0.4 
0.05±0.01 
8.0 
1.2/45 


4.6±0.9 
1.26±0.06 

18±2 
1.07±0.04 
6.8±0.2 

0.6/ 
180±60 

1.0 

0.8/47 


3.8±0.8 
1.21±0.06 

18±2 
1.34±0.08 
6.6±0.1 

0.6' 
160±50 

1.5 
0.6/47 


2±1 
0.7±0.1 

15±2 
5.4±0.9 
6.4±0.2 

0.4' 
110±60 

18.2 
1.1/47 




INTEGRAL/ISGK1 + RXTE/PCA energy spectra* 






Photon index 

E mt (keV) 

PL norm* 

Flux c (2-100 keV) 

X 2 r /dof 




0.81' 
21±1 
1.5±0.9 
10.7 

1.1/63 




1.21' 
19±4 
0.2±0.2 
1.9 

0.7/62 


0.7' 
22±1 
0.7±0.3 
27.1 
1.1/62 


Power spectra 


Power-law index 
rms (%) (0.01-1 Hz) 


1.7+0.1 
24±3 


1.6±0.1 

25 ±3 


2.1±0.4 
18±5 


2.2±0.4 
18±4 


1.39±0.05 
25+1 



**: It includes the model to the RXTE/PCA spectra with the same value of the parameters 
a: xlO 22 cnT 2 c: xl0~ 10 erg era" 2 s~' 

b: xl0~ 2 cm" 2 s _1 keV" 1 at 1 keV /: fixed 



this new component to the best-fit value found in the fit to the PCA 
data and allow only the cutoff energy and normalisation to vary. 
In the faint state the addition of INTEGRAL data to the PCA fit 
represents a continuation of the lower energy spectrum, that is, the 
two cutoff energies are consistent within the errors. In contrast, the 
broad-band spectra for the bright and flare state are not well fit- 
ted by a single power law plus cutoff energy. The extension of the 
spectra above 30 keV for these two states requires an extra cutoff 
energy. Table [4] shows the results of the spectral fits. The parame- 
ters that fit the PCA spectra also fit the PCA+INTEGRAL spectra. 
The only difference is the extra exponential decay component. 

The large effective area of the RXTE instruments make it espe- 
cially valuable for timing analysis of intensity variations from high- 
energy sources. Therefore the study of the aperiodic variability of 
SAX J2 103. 5 +4545 was performed with PCA data only. Given the 
large number of observations and in order to obtain an homoge- 
neous set of power spectra in terms of energy and frequency range 
coverage we used the configuration mode Standard 1. The light 
curve corresponding to each PCA observation was binned with bin 
size 0.125 s, which is the maximum resolution provided by the 
Standard 1 data mode, and divided into segments of 128 s. An FFT 
was calculated for each segment giving power spectra covering the 
frequency range 1/128-4 Hz. The final power spectrum for each 
observing interval resulted after averaging all the individual power 
spectra and logarithmically rebinning in frequency (Fig.[6jl. In gen- 
eral, the power spectra were well fitted to a single power law. The 
values of the power-law index and fractional amplitude of variabil- 



ity (rms) given in Table|4]correspond to the mean and standard de- 
viation of the average values from all the intervals included in each 
state. Given the intrinsic variability of the source, we preferred to 
obtain mean values from the relevant observations rather than one 
single average power spectrum from the entire light curve of each 
state . Intervals resulting with less than six segments (i.e. six power 
spectra) were not considered because of poor statistics (low num- 
ber of points in the low-frequency bins). Likewise, power spectra 
whose fits to a power law gave reduced x 2 larger than 2 were re- 
jected. Since the flare state include one only observational interval 
the values quoted are those from the average power spectrum of that 
interval. From the results shown in Table|4]we conclude that the red 
noise component becomes flatter as the X-ray flux increases. 



3.4 X-ray spectral-temporal correlation in the bright state 

In addition to characterising the bright and faint states we also 
searched for correlations between the temporal and spectral pa- 
rameters on shorter timescales. For each PCA observation falling 
in one of the two bright or two faint states defined above we ob- 
tained an energy spectrum and two power spectra corresponding to 
the energy ranges 3-10 keV, 10-30 keV, respectively. In this case 
the light curve at different energies were extracted from data of 
the GoodXenon configuration mode. As before, the 3-30 keV en- 
ergy spectra were fitted with an absorbed power law plus exponen- 
tial cutoff and a Gaussian representing fluorescence from cold iron, 
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Figure 5. Joint PCA+ISGRI spectra for the flare, bright2 and faint2 states. 
The model components comprises two cutoff power laws, a Gaussian and 
low-energy absorption. Circles correspond to PCA data and squares to IS- 
GRI data. 
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Figure 6. Power spectra for the flare (MJD 54221), bright2 (MJD 52653) 
and faint2 (MJD 53265) states. 
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Figure 7. Relationship between the power-law index and rms resulting from 
the fits to the power spectra in two energy bands for the two bright states. 
Circles represent state 1 and squares correspond to state 2. 
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Figure 8. The Tspec - Tpsd diagram. Slope of the power spectral continuum 
as a function of the photon index of the energy spectrum. Circles represent 
state 1 and squares state 2. For clarity the error bars of the faint2 state were 
omitted. 
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Figure 9. Photon index and power-law index as a function of X-ray flux in 
the 3-30 keV range. Circles represent state 1, while squares correspond to 
state 2. For clarity the error bars of the faint2 state were omitted. 



while the power spectra exhibited strong red noise that was fitted 
with a single power law. 

We first checked whether the aperiodic variability, i.e., the 
continuum of the power spectra, depended on energy. Fig.|7]shows 
the power-law index obtained from the 3-10 keV and 10-30 keV 
power spectra for the bright states. A fit of both sets of data to a 
straight line gives a slope of 0.83+0.06, indicating that the power 
spectral continuum weakly depends on energy. That is, on aver- 
age, the slope of the power spectra at higher energies is flatter than 
at lower energies. The fractional amplitude of variability in the fre- 
quency range 0.01-1 Hz is also, on average, higher at high energies. 
As it can be seen in Fig. [7] the rms data points measured from the 
10-30 keV power spectrum lie systematically above the one-to-one 
relationship, denoted by the straight line. 

Then we investigated the interdependency of the spectral and 
temporal parameters. Figure [8] shows the relationship between the 
power-law index resulting from the power spectra using the 2-60 
keV (Standard!) light curves (Fpso) and the photon index from the 
energy spectra (Fspec)- It is clear that in the bright state steeper 
energy spectra correspond to steeper power spectra. However, this 
correlation is lost in the faint state. 

The variation of the temporal and spectral parameters as a 
function of the 3-30 keV flux is shown in Fig. [9] Again cor- 
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4 DISCUSSION 



4.1 Optical analysis 



0.5 

Orbital phase (P 



= 12.66528 d) 



Figure 10. Spectral and timing parameters as a function of orbital phase for 
the bright state. Circles correspond to state 1 and squares to state 2. 



relations between these parameters are apparent in the bright 
state only. The TspEc-flux plot confirms the known result that 
the source spectrum becomes har der for higher flux (Bav kal et al.l 
l2002l :l Camero-Arranz e t al. 2007). However, we find that this cor- 
relation also holds in the timing domain with flatter power spectra 
being associated with higher flux. The r SPEC - r PS D (or r PSD -flux) 
represent new results, not previously reported on a high-mass X-ray 
binary. 



3.5 Orbital phase variability 

We also searched for orbital variability in the bright and faint states. 
We assigned an orbital phase to each PCA observation falling 
within the four intervals defined in Table [4] using the orbital solu- 
tion of ICamero-Arranz et al.l J2007h (no signi ficant difference was 
found in the results if the orbital parameters of lBavkal et alj fc007h 
were used). The data points were then grouped into five phase bins. 
Figure [TOl shows the results of the orbital variability analysis. We 
show the bright state only as no correlation with orbital phase was 
found in the faint state. The X-ray flux in the 3-30 keV band in units 
of 10~ 10 erg cnT 2 s , the photon index of the energy spectra, Tspec, 
the power-law index of the power spectra, r PSD , the fractional am- 
plitude of variability in the frequency range 0.01-1 Hz, rms, and 
the intensity of the iron line, /p c , are shown as a function of orbital 
phase. In this figure, periastron corresponds to orbital phase ~0.4, 
while maximum flux is reached around phase 0.7. 

When the 3-30 keV flux is high, the intensity of the iron line is 
also high, while the photon index and power-law index are low. The 
rms also shows a weak modulation: the source appears more vari- 
able when the flux is high, but the relatively large errors diminish 
the significance of this trend. The iron line energy remains fairly 
constant throughout the orbit. A weighted mean gives 6.44+0.06 
keV. 



Numerous studies on Be stars toachs et all 1 1984 



Hummel & Vranckenl [l995l: iHanuschikl 1 19961 : I Porter & Riviniuj 



120031 : iTvcner et al] 120051 ; Ijones et al.l |2008|) showed that the 
equatorial disc in Be stars is rotationally supported consistent 
with a Keplerian velocity field and that the strength and shape 
of the Ho- line provides important information on the structure 
of the disc. Double-peaked profiles have been interpreted as a 
result of the Keplerian motion of the particles in the disc, with the 
peak separation giving an estimate of the rotation velocity. Split 
profiles are expected to occur when the disc is smaller, that is, at 
smaller radii. When the disc grows, since the Keplerian velocity 
is inversely proportional to the radius, the velocity becomes too 
small to be discernible in low or intermediate resolution spectra. 
This trend of double-peaked profiles merging in to one single-peak 
line as the Ho- equivalent width i ncreases is a well known re sult 
both in Be stars foachs et all 19861) and Be/X teeig et alliooch . 

We find a good correlation between the strength and shape of 
the Ho- line and the X-ray intensity in SAX J2103.5+4545 (Fig.©. 
Emission line profiles are seen only when the source is bright in X- 
rays, while low-intensity X-ray states are associated with absorp- 
tion line profiles. This type of correlation has been seen in other 
Be/X. What makes SAX J2103.5+4545 unique is the rapidity of the 
changes in the strength and shape of the Ha line, presumably due 
to a the rapid growth of the disc and the almost immediate trans- 
fer of matter into the neutron star to generate X rays. Typically, the 
time scales for the formation and disappearance of the circumstel - 
lar disc in BeX are larger than 3 years (see e.g. iReig et alj .2005). 
In SAX J2103. 5+4545 it is 1-2 years. Thanks to our monitoring 
in the optical band we can put some constraints on the duration of 
the formation and subsequent loss of the disc (see Fig. [T] and Ta- 
ble[TJ. The Ha line was in absorption with EW(Ha)=+2.4 A with- 
out any trace of fill-in emission in June 2006. The largest value of 
the EW(Ha), which would roughly correspond to the maximum ex- 
tension of the circumstellar disc, was achieved in September 2007, 
that is, SAX J2103.5+4545 took about 14 months to develop the 
circumstellar disk. This duration should be considered as an up- 
per limit. Given the drop in X-ray intensity around January 2007 
and the lack of observations in the optical band at this time, it is 
very likely that the Ha was in absorption at the beginning of 2007. 
If this is the case then the formation of the disk took just about 
8 months. By December 2007, the strength of the Ha line was 
clearly in the decline, turning into an absorption phase (EW(Ha) 
had weakened to just 1 A, although still in emission). Full absorp- 
tion at a level of EW(Ha)=+2.3 A was seen about one year af- 
ter maximum. We conclude that the total duration of the forma- 
tion/dissipation of the circumstellar disc in SAX J2103. 5+4545 is 
< 2 yr, and most probably of the orde r of 1.3-1.5 yr. T his dura- 
tion compar es to 3-5 yr of 4U 1 15+63 dReig et alJioOTT). 4-5 vr o f 
A 0535+26 taaigh et al.ll2004l) a nd V0332+53 jGoranskiill200ll) 
4 yr of LS992 JReig et al.ll200ll) and 7 yr of X-Per faark et al.l 
2001). The reason for these fast changes must lie in the relatively 
narrow orbit (SAX J2103.5+4545 is the Be/X with the shortest or- 
bital period). 

In SAX J2103. 5+4545 changes in the disc happen so fast that 
the system does not have the time to develop a large disc. The fact 
that i) the largest EW(Ha) measured in SAX J2103.5+4545 is ~ 
-5 A, a relatively small value and ii) the emission line profile ap- 
pears always as double peaked, provide evidence in support of a 
small disc radius. The narrow orbit and continuous passages of the 
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neutron star limit the growth of the disc and do not allow the Be 
star to develop an extended disc. On the other hand, the duration 
of the disc-loss phase also seems to be short. In fact, it is likely 
that the complete loss of the disc never materialises. An inspection 
of Table Q] reveals that CQE and emission profiles are present in 
the majority of the observations. As stated above emission lines are 
formed in the disc. CQE profiles also require that a circumstellar 
disc e xists, although of a much smaller spatial extent dRivinus et al.l 
1999). Note also that some of the profiles classed as pure absorp- 
tion may contain a weak central peak that is not apparent due to 
the low spectral resolution. Another interesting characteristic of the 
X-ray behaviour of SAX J2103. 5+4545 is the large and sudden in- 
crease in intensity that precedes some of the outbursts. The X-ray 
intensity increases by a factor of ~20 in a few days. The onset of 
these flares occurs at or near periastron and last for no more than 
two or three orbital cycles. Then the X-ray intensity abruptly falls 
to almost pre-flare values. Immediately after, a lower intensity and 
longer outburst begins (the bright state). 

Thus the picture that emerges to explain the long-term X- 
ray/optical variability is the following: the Be star's equatorial disc 
forms due to ejections of matter from the Be star's photosphere 
through a still unknown mechanism. As the disc grows the Ho- line 
turns from absorption into emission. The fact that the emission line 
profile is always split, i.e., no single-peak line is seen and the rel- 
atively small EW(Ha') indicates a small disc. We can estimate the 
disc radius from the separation of the Ho- profile peaks, as this value 
and the outer radiu s (Rdisc) of the emission line forming region are 
related by (see e.g. lHummel & Vrancken|[l995t) 

^disc /2vsinA~ 
R » \ A pcak / 

where v sin i is the projected rotational velocity of the B star (v is 
the equatorial rotational velocity and i the inclination toward the 
obser ver). For SAX J2103.5+4545 v sin i » 240 km s" 1 dReig et al.l 
l2004h . while the average peak separation, as measured from the 
September 2007 spectra is A peak = 340 km s _1 . Thus R disc ~ 2 - 3 
R*. Taking the characteristic values of mass and radius of a B0 star, 
name ly M, = 20 M Q and R, = 8 Rg for the B star IVacca et all 
1996) and the canonical mass of a neutron star M x =1.4 M Q , the 
orbital period (P OI b = 12.6 d) implies, according to Kepler's third 
law, an orbital separation of a a 63 Rq or equivalently a as 8 R,. 
Correspondingly, periastron lies at R per = 4.7 R, for an eccentricity 
e = 0.4. Similarly, the radius of the Roche lobe of the primary star 
is dPaczvnskil 197 ll) R„, la = 0.38+0.20 log(M»/M.,) = 0.6. That is, 
the radius of the disc is similar to the size of the Roche lobe, which 
in turn, is similar to the periastron distance, R RL w R pa ~ /?disc- 

We conclude that due to the small orbit, the neutron star trun- 
cates the disc at very small radii. It is very likely that, occasion- 
ally, the neutron star physically impinges on the disc during peri- 
astron passage giving rise to sharp and intense X-ray flares. Note 
that the fact that R iisc = 2R, does not mean that the disc terminates 
abruptly at that distance, but matter would extend to larger radii 
following a density law of the form p ~ r~P , where fj w 2.5 - 3 
(Watersj fl986l Il988l) . The fact that the optical indicators (magni- 
tudes and EW(Ha)) continue rising after the flare episode implies 
that the disc has not yet been destroyed. In subsequent periastron 
passages, matter is still transferred to the neutron star but in a less 
dramatic way. An accretion disc around the neutr on star is likely to 
be formed, as indicated by the spin-up episodes l ilnam et alii 20041 ; 
ICamero-Arranz et al.l l2007). When the fuel that power the X-ray 
activity is used up, that is to say, when the circumstellar disc dis- 



appears then accretion ceases and the source returns to an X-ray 
quiescent state. 

4.2 X-ray analysis 

Our X-ray spectral analysis of the bright a nd faint states of SAX 
J2103.5+4545 confirms pre vious reports dBavkal et alj200ll2007l : 
ICamero-Arranz et al.ll2007f) that in the bright state the spectrum is 
harder than in the faint state. This result is opposite to what it is 
seen in other types of X-ray binaries but it is not exclusive of SAX 
J2103. 5+4545. Black-hole syste ms and low-mass X-ray binaries 
display different source states dvan der Klis 2006, and references 
therein). In black-hole systems three basic states are identified, 
namely, soft, hard and intermediate (which includes the so called 
very high state). Although these states may occur at any luminos- 
ity, when the black-hole system is followed through an outburst, 
the soft state normally corresponds to higher X-ray intensity val- 
ues, while the hard state is normally seen at the beginning and end 
of the outburst, i.e, at lower count rates. Likewise, spectrally hard 
states in low-mass X-ray binaries occur at lower luminosities tha n 
softer spectra states dOlive et alJl2003l ; lMaccarone & Coppill2003l) . 
In contrast, in a spectr a l and timing study of four Be/X-ray during 
major outbursts, iReid ( 2008) found that three of the sources in- 
vestigated showed low/soft spectral states (the so-called horizontal 
branch) at the beginning and at the end of the outburst. 

The origin of these differences must lie on the accretion pro- 
cesses in the vicinity of the compact object. Unlike black-hole sys- 
tems and low-mass X-ray binaries with weak magnetic-field neu- 
tron stars, the compact object in Be/X-ray binaries is a strongly 
magnetised neutron star, which results in signi ficant coupl ing of 
the magnetic field and the accretion flow l l White et alj [l983j). One 
of the results of this coupling is that the accretion flow is threaded 
onto the neutron star magnetic field lines and channeled on to 
the magnetic poles, producing two or more localised X-ray hot 
spots (see e.g iFrank et alj|2002l) . Above these hot spots the mag- 
netic field lines adopt the configuration of a roughly column- 
shaped surface or funnel, which is referred to as the accretion 
column, the polar cap being its base. The polar caps are treated 
as dense thermal mounds with a blackbody spectrum. The detec- 
tion of such blackbody components at low energies in a number 
of low-luminosity X- ray pulsars give support to the existence of 
such thermal mounds dCoburn et al. 2001 ; Mukhe riee & Paull2005l : 
lLaPalombara&Mereghettill2006l ; lReigetalJl2009r) . The surface 
of the mound is respo nsible for photon creation and absorption 
dBecker &~W olff 2005). The high-energy emission observed from 
accreting pulsars results from the formation of a standing shock 
wave on this accretion column: the energetic particles in the accre- 
tion fl ow comptonise the t hermal emission produced on the polar 
cap dBecker & Woiffll 20051) . Increased mass accretion rates are ex- 
pected to result in hard er X-ray spectra dLanger & Rappapor J 19821 : 
Becker & Wolff 20051) because the photons spend more time, on av- 
erage, being upscattered in the flow before escaping. 

The faint state corresponds to a quiescence and stable state 
of the source in which the X-rays are no longer powered by disc 
accretion but most likely are the result of wind accretion. Wind- 
fed system, like supergiant X-ray binaries, show erratic and flaring 
(high-amplitude changes on timescales of seconds) X-ray variabil- 
ity that it is consequence of inhomogeneities in the accretion wind. 
This type of variability is almost absent in SAX J2103. 5+4545, 
where the main source of variability at short timescales are the X- 
ray pulsations. The fact that we do not see random high-amplitude 
variations during the faint state of SAX J2103. 5+4545 might in- 
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dicate that the dominant accretion wind is the more stable equato- 
rial lo w-velocity high-density wind characteristic of BeX {Waters 
1988). This wind, that ultimately would create the circumstellar 
disc, would be highly suppressed during the faint state. It would 
also indicate that the neutron star orbit is coplanar to the plane of 
the circumstellar disc. Otherwise, the neutron star would be ex- 
posed to the higher velocity and more inhomogeneous polar winds, 
typical of massive stars, giving rise to erratic X-ray variability. The 
narrow orbit of SAX J2103. 5+4545 would also help avoid the po- 
lar winds if the orbit is coplanar to the disc. Occasionally sudden 
and short-lived increases in X-ray flux are seen superimposed on a 
long-term decline phase during the faint state (see e.g. MJD 52800- 
53200 in Fig.gJ. These mini-outbursts may correspond to fail at- 
tempts to create the circumstellar disc, i.e. ejection of matter from 
the Be star's photosphere that are not retained by the star. 

In contrast, the bright state exhibits high richness in variability. 
Both spectral and temporal parameters correlate with one another 
and with X-ray flux and orbital phase. During the orbital maximum 
the spectral photon index is minimum. Similar variability behaviour 
has been reported in the supergiant X-ray binary 2S 0114+650 
dFarrelletal.ll2008h . Based o n the similarities with the eclipsing 
X-ray pulsar EXO 1722-363 ^Thompson et alJl2007l) . lFarrell et all 
( 2008) attributed the spectral variability to absorption effects as the 
neutron star passes behind a heavily absorbing wind. 

We have found for the first time in a HMXB a correlation be- 
tween the shape of the energy and power spectra, namely softer en- 
ergy spectra correspond to the disappearance of power at high fre- 
quencies, i.e. softer power spectra. Although there are no detailed 
models of the power spectra arising from accretion onto highly 
magnetized neutron stars, one would expect that the simple sce- 
nario described above would also explain the Tspec - Tpsd corre- 
lation (Fig. [8} At high mass accretion rates, i.e., when the X-ray 
flux is high, the magneto sphere shrinks, as the m agnetosphere's ra- 
dius r m oc Mr 111 (see e.g. lDavidson & O striker 1973j). If the source 
of variability in the power spectrum comes from the interaction 
between the accretion disc and the magnetosphere then smaller 
magnetospheric radius implies shorter characteristic time scales (or 
higher frequencies), hence flatter power spectra. When the flux de- 
creases (softer energy spectrum), the size of the magnetosphere in- 
creases and so does the characteristic timescale of variability pro- 
ducing more power at low frequencies, hence steeper power spec- 
tra. The fact that the rms is similar in the bright and faint states 
supports the notion that the source of variability is similar in both 
states and the only difference is the size of the emitting region. 

If we assume that the characteristic frequencies result from 
Keplerian motion in the inner parts of the accretion disc, and that 
the accretion disc extends all the way down to the magnetosphere, 
then the characteristic Keplerian frequency of an orbiting free par- 
ticle is 



GM 

Assuming that all the kinetic energy of infalling matter is given 
up to radiation, i.e., L x = GM X M/R X , an increase of the X-ray 
luminosity by a factor of ~6, as observed in the bright state of 
SAX J2103. 5+4545, implies a reduction in the size of the magneto- 
sphere by ~ 40%, while the characteristic frequency of variability 
increases by a factor ~2. 

In Fig.|7]we showed that the fractional rms is larger at higher 
energy. If the source of variability were variations in the soft pho- 
ton input, then simple Comptonization models would predict a de- 
crease of the amplitude of variability with increasing photon en- 



ergy. This trend is the result of the averaging effect of Comp- 
ton scatterings: low-energy photons retain most of the variabil- 
ity of the seed photon input because they have not suffered many 
scatterings. In contrast, the variability of high-energy photons is 
smeared out as these photons spend longer time in the comptoniz- 
ing medium. These prediction are not consistent with observa- 
tions of black-hole systems. In low-mass X-ray binaries and black- 
hole systems, the fractional rms of discrete noi se components 
(i.e. QPO) tend to increa se with photon energy (Gilfan ov et al.1 
120031: lHoman et al.1 l200l1) . In the case of broad-band noise, the 
relationship of the rms with e nergy shows a much more com- 
plex behaviour Jlin et alJb oOO). and it is source state dependent 
( IGierlinski & Zdziarskil 120051) . Note also that in black-hole bina- 
ries the geometry of the Comptonizing medium is very different 
from the geometry in the pulsar accretion column and the optical 
depth is expected to be lower. In order to have larger fractional rms 
at higher energies, as shown in Fig[7] the source of variability can- 
not be solely due to intrinsic variat ions of the source of soft pho- 
tons. Gieiiinski & Zdziarski l l2005t) found that an increase in the 
fractional rms with photon energy can be achieved if the power re- 
leased in the comptonising plasma changes, i.e., if the energy of the 
comptonising electrons changes (electron acceleration and/or direct 
heating). Applied to SAX J2103.5+4545 this would mean that the 
power supplied to the electrons in the accretion flow changes, per- 
haps as a result of varying mass accretion rate, which may well 
affect the geometry and optical depth of the Comptonizing column. 



5 CONCLUSIONS 

We have performed an analysis of the X-ray timing and spectral 
properties of the Be/X-ray binary SAX J2103. 5+4545 in correla- 
tion with optical spectroscopic observations of the Ho- line. We find 
a good correlation between the strength and profile of this line and 
the X-ray activity of the source. Emission line profiles are observed 
during bright X-ray states. In these states the Ho- equivalent width 
is largest. During faint X-ray states Ha appears in absorption. SAX 
J2103. 5+4545 exhibits the fastest time scales for the disappearance 
and reformation of the circumstellar disc, compared to other typi- 
cal Be/X-ray binaries. In less than 2 years SAX J2103. 5+4545 is 
capable to form and lose the disc. Due to its narrow orbit, the neu- 
tron star prevents the disc from extending to large radii; the disc 
is truncated at 2-3 stellar radii. Based on the Ha line profile we 
set stringent constraints on the size of the circumstellar disc and 
we find that its size can reach the Roche-lobe radius of the system 
(which it is similar to the periastron distance of the two objects). 

We have also performed for the first time in a Be/X-ray bi- 
nary a correlated study of the X-ray spectral and timing properties 
and found a correlation between the shape of the energy spectrum 
and that of the power spectrum in the bright state. When the energy 
spectrum is soft, the power at high frequencies is suppressed. This 
correlation, together with the existence of hard/high and soft/low 
states, has been interpreted in the context of emission from the ac- 
cretion column. 
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